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been computed with the help of (2) - obtained from (L, L r ) = (16,10), (16,6) and (10,6) are 
displayed in Tablet. Again this may be compared to 0.67 and 0.77 from first-order and 
second-order e-expansion. 



0.1 

m(t) 

0.09 
0.08 
0.07 
0.06 




50 100 150 200 



t (MCS) 

Fig. 4: Test of finite-size scaling relation (3). Compared are profiles for L = 16 and 
mo = 0.06 (solid line) and L = 10 and mo = 0.082, 0.084, 0.086 (dashed lines, from bottom 
to top). For details see text. 

In summary, we have investigated the influence of the initial conditions on the relax- 
ational behavior of the three-dimensional Ising model by means of Monte Carlo simula- 
tion. The measured values of the exponents 8' and xq are consistent with those found by 
renormalization-group improved perturbation theory [1]. As demonstrated in Fig. 3, also 
the dependence of the amplitude of the linear decay on the initial magnetization agrees 
very well with the expectation from continuum calculations [4, 5] 

Compared with what is standard today in Monte Carlo simulations in critical dynamics 
[8, 9], the lattices we have studied are relatively small. However, we think the results 
derived from our small systems are very well suited as a first step to verify the anomalous 
initial behavior and other universal properties related to it. Especially, we point out that 
the dependence of 6' on the lattice size in the range of L studied is very weak (cf. Table 
1). Thus, we do not expect qualitative or dramatic quantitative changes for the simulation 
on larger lattices. 



8 



of Fig. 2. The results which are displayed in Fig. 3 clearly show the linear increase for small 
mo and the tendency to become independent of mo for larger values. Fitted to the numer- 
ical data is the analytic large-n solution (9) with constants A = 1.98 and B = 14.2. We 
actually do not expect that the large-n result also holds in exactly the same form for the 
Ising system, but apparently the scaling function for the latter is quite similar. 
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Fig. 3: Amplitude of the linear decay as a function of mo for L = 10. The squares resp- 
resent the Monte Carlo results. The function fitted to them is the exact large-n solution (9). 

Finally, we test the scaling form (3) and determine directly the scaling dimension xq by 
comparing profiles from systems of different size. If (3) is satisfied, it should be possible to 
map curves onto each other by appropriate rescalings of mo, m and t. To this end, we have 
calculated one profile for lattice size L and initial magnetization mo and several curves for 
a smaller lattice L' and a certain range of mj. Now the curves from the primed system 
can be mapped to the unprimed system by (i) multiplying by an overall factor b^l v and (ii) 
by stretching the time axis by b z with b = L/L' . Here we have used the literature values 
z = 2.04 [8] and (3/v = 0.52 [11], which are, by the way, consistent with our own data. For 
one example the result is shown in Fig. 4. The solid line is from L = 16 and mo = 0.06, and 
the three dashed lines are those from L' = 10 and m'g = 0.082, 0.084 and 0.086, respectively. 
(Note that due to the finite number of the degrees of freedom the increment for changes of 
m Q here is 0.002.) After matching the magnetizations by linear interpolation of m Q , xq is 
determined with the help of m'g = (L/L') x ° rriQ. The results for xq and 8' - the latter has 
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Fig. 2: Magnetization profiles for L = 10 and mo = 0.02, 0.06, . . . , 0.38 in semi-logarithmic 
representation. 

Because we expect the universal initial increase only for small mo, we calculate 8' from 
eight profiles with 0.02 < mo < 0.10 and 4 < L < 16. The results and their statistical 
errors are displayed in Table 1. (The figure for mo = 0.04 and L = 10 is included in the 
table but the corresponding curve is not shown in Fig. 2.) They have to be compared with 
0.08 from first-order and 0.13 from second-order e-expansion [1]. In the range of L and 
mo we consider, the scattering of the data is small. Closest to the truth is probably the 
6' = 0.102(2) from L = 16 and m = 0.06. 
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8' 


0.02 


0.115(11) 


16 


0.102(2) 


(16, 6) 


0.725(13) 


0.100(6) 


0.04 


0.105(5) 


12 


0.102(3) 


(10, 6) 


0.738(26) 


0.107(13) 


0.06 


0.106(5) 


8 


0.114(10) 


(16, 10) 


0.710(22) 


0.093(11) 


0.10 


0.102(6) 


4 


0.096(43) 





Table 1: Monte Carlo results for 9' and xq. 



The dependence of the amplitude of the linear decay has been obtained from the data 



6 




Si = — 1. Afterwards, if the avarage magnetization is below or above the desired value mo, 
single spins are chosen randomly and flipped if Si = — 1 or Si = +1, respectively, until the 
average magnetization is exactly rriQ. 

Then the standard heat-bath formalism is used to trace the time evolution of the system 
[7]. Time is measured as usual in Monte Carlo steps per spin (MCS) and measurements 
are carried out after each sweep. We work at the (bulk) critical temperature J/T c = 0.221G 
[11]. For each magnetization profile we have averaged over 150,000 to 400,000 statistically 
independent histories of the system or runs, up to t = 100 — 200 MCS. Further, for each 
run a new initial configuration is generated and a new sequence of random numbers is used 
for the ensuing relaxational process. 

Profiles for various L and mo are displayed in Fig. 1 and 2. In Fig. 1 both axis are log- 
arithmic to show clearly the power-law behavior (5) for early times, whereas in Fig. 2 only 
the m-axis is plotted logarithmically to bring out the linear decay for long times. All curves 
show the initial increase, and 6' is determined by fitting a straight line to the short-time 
regime. 
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Fig. 1: Magnetization profiles for mo = 0.06 and L = 4, 8, 12, 16 in double-logarithmic 
representation. 
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Ritschel [4]. For the special case of fixed L, their result for the amplitude takes the simple 
form 



with non-universal constants A and B. For the Ising system (n = 1) we expect a qualita- 
tively similar behavior of the asymptotic amplitude, i.e., the linear dependence for small 
mo and asymptotically independence of the initial magnetization for large mo. 

Now, the question asked in this Letter is: Do the results on universal short-time behavior 
reported above stand up to the Monte Carlo "experiment"? 

Much numerical work has been devoted to the simulation of time-displaced correlation 
functions in equilibrium [8, 9], but to our knowledge the process described above has not 
been studied in the literature so far. However, projects similar to ours are at present under 
investigation [10]. For the time being, our major aims may be summarized as follows: 

• Simulation of the short-time behavior of the relaxational process starting from the 
non-equilibrium initial state as explained above and measurement of the exponent 6' . 

• Measurement of the universal amplitude of the exponential decay for late times and 
comparison with (9). 

• Test of the scaling form (3) and thereby determine the scaling dimension xq. 

An initial state with non- vanishing magnetization and small correlation length in princi- 
ple can be generated by letting evolve the system to thermal equilibrium with temperature 
T T c within an external magnetic field. However, this procedure turns out to be unprac- 
tical because it is too time consuming, and, more importantly, due to the finite number 
of degrees of freedom, even for T > T c and strong magnetic field, fluctuations of the ini- 
tial magnetization of the canonical ensemble can not be neglected. As a consequence, the 
width of the distribution of the initial magnetization gives rise to an additional scale which 
may alter the simple relation (3). (In the continuous system the same phenomenon occurs 
if the initial temperature is close to T c [1, 4].) In order to avoid this, we start from a 
"microcanonical" initial state with fixed mo and small correlation length. An initial con- 
figuration with a given mo is generated as follows. First, the spins on single lattice sites 
are determined independently, with probability 1/2 + mo/2 for Si = 1 and 1/2 — mo/2 for 
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(9) 
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where the initial time scale to is given in (1), and ti is the well-known finite-size relaxation 
time ti ~ L z [6]. 

Thus, as stated in (4), the relaxational behavior of the critical Ising system starting from 
a non-equilibrium initial state with small correlation length (<L and initial magnetization 
mo is governed by two macroscopic time scales. This gives rise to clearly distinguishable 
stages of the relaxational process. First, the characteristic initial increase of order, 

m(t) ~ t 6 ' , (5) 

discovered in [1] for the infinite system, is also found in the finite system [4] when t < to, tj^. 
The only condition for the increase is that to remains macroscopic which, in turn, means 
that mo has to be small (c/. (1))- Second, at time t max the magnetization reaches its 
maximum and then starts to decrease towards the equilibrium value m(t —t oo) = 0. For 
ti S> to, the time t max is roughly given by to, but t max decreases with decreasing ti and can 
be much smaller than to for t^ -C to- Third, the long-time behavior for t ^> t^ is described 
by an exponential 

m(t > t L ) a Moo e _i/ti . (6) 

While the time scale in this "linear" regime is not influenced by the initial state [5], the 
amplitude M m is a universal function of both mo and L. It may be written in the form [5] 

Moo(L, m ) « const. L~^ v Q{to/t L ) . (7) 

In the limit tj, ^> to, the scaling function Q approaches a constant. Hence, in this case 
the dependence on to (and therefore on mo) drops out of the long-time behavior of m(t). 
On the other hand, for t^ -C to the scaling function behaves as Q(w) a w~ x °/ z . Thus, the 
amplitude becomes [4] 

Moo(L, mo) « const, mo L z6 ' , (8) 

i.e., it is proportional to the initial magnetization mo, and the exponent 8' that in the 
infinitely extended system only appears in the short-time regime now governs the in- 
dependence of Mqo. In other words: In the finite system there is a long-time memory 
of the initial condition in the sense that the asymptotically leading term depends on it. 

For the Ising model the scaling function Q has not been calculated so far. However, for 
the general n-vector model an exact solution for n — > oo has been obtained by Diehl and 
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Five years ago Janssen et al. [1] discovered that under certain circumstances relaxation 
processes in Ising-like systems near criticality show strong anomalies. When for example a 
system belonging to the dynamic universality class of model A [2] is quenched from an state 
with temperature T > T c to the critical point, an initially existing small magnetization will 
grow for a macroscopic time span before the actually expected decay towards equilibrium 
takes over. The initial increase of order, which was later found also in other dynamic models 
[3], is described by a universal power law with a new dynamic exponent, 9' , that can not 
be expressed in terms of a scaling relation between the static exponents and the dynamic 
exponent z. The reason for the short-time anomaly has to be sought in the initial conditions, 
more precisely, in the behavior of the initial magnetization mo under renormalization-group 
transformations. As found by Janssen et al., the operator mo has a scaling dimension xq 
that is in general different from the dimension of the equilibrium magnetization, = j3 jv. 
As a detailed scaling analysis reveals [4], this gives rise to a macroscopic time scale 

k ~ m - 2/ *° , (1) 

and the exponent 8' is related to xq by 

6>= X -^±. (2) 

z 

Numerical values for the new exponents have been obtained by means of the e-expansion 
[!]• 

More recently, Diehl and Ritschel [4] extended the above-mentioned considerations to 
systems of finite size. In the following these results which have been obtained by finite size 
scaling analysis and with the help of continuum field theory will be summarized because 
they are essential for the interpretation of the Monte Carlo data. The basic quantities that 
enter the scaling analysis are the time t, the lattice size L in units of the lattice constant 
a, and the spatially averaged magnetizations mo and m(t). At bulk criticality and for L 
and t much larger than any microscopic scales, one expects the following scaling form of 
the magnetization: 

m(t, L, m ) ~ l~ 13/11 m(r z t, l^L, l x ° m ) . (3) 
By choosing the arbitrary rescaling parameter / ~ t 1/z , one can rewrite (3) as 

m(t, L, m ) ~ t-M" z F(t/t L , t/t Q ) , (4) 
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Abstract: The time evolution of the three-dimensional critical Ising model relaxing from 

a non-equilibrium initial state is studied by means of Monte Carlo simulation. We observe 

g the characteristic initial increase of the (spatially) averaged magnetization predicted by 
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Janssen et al. [1]. The exponent 8' that governs the initial behavior is determined, and the 
O dependence of the long-time linear decay on the initial magnetization, mo, is analyzed. Our 

simulation corroborates earlier results derived from continuum models. 
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